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Abstract 

We have derived a new set of semiclassical equations for electrons in mag- 
netic Bloch bands. The velocity and energy of magnetic Bloch electrons are 
found to be modified by the Berry phase and magnetization. This semiclas- 
sical approach is used to study general electron transport in a DC or AC 
electric field. We also find a close connection between the cyclotron orbits in 
magnetic Bloch bands and the energy subbands in the Hofstadter spectrum. 
Based on this formalism, the pattern of band splitting, the distribution of 
Hall conductivities, and the positions of energy subbands in the Hofstadter 
spectrum can be understood in a simple and unified picture. 

PACS numbers: 72.10.Bg 72.15.Gd 72.20.Mg 
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I. INTRODUCTION 



The semiclassica method has played a very important role in studying electron dynamics 
in periodic systems.0 In this approach, the effect of a periodic potential is treated by quantum 
mechanical methods and yields usual band structure for energy spectrum, while an extra 
electromagnetic field is treated as a classical perturbation. The velocity of an electron in 
the one-band approximation is given by 

where £n is the energy spectrum for the n-th band. The dynamics of quasi-momentum k is 
governed by the Lorentz-force formula 

/ik = -eE - er X B, (1.2) 

where E and B are the external electric and magnetic fields. These equations may be 
regarded as the equations of motion for the center of mass of a wave packet in the r and k- 
spaces. Tremendous amoimt of work has been done to justify these simple looking formulae 
and to their quantization.B 

These formulae, however, become invalid if the magnetic field is so strong that it is 
no longer appropriate to be treated as a perturbation. In this case, we need to solve the 
Schrodinger equation with the following Hamiltonian: 

i^o = ^(-^/^|: + eAo(r)) + ^r), (1-3) 

where Ao(r) is the vector potential of a homogeneous magnetic field,il and V{y) is a periodic 
potential. The eigen-energies of Eq. (1.3) will be called magnetic Bloch bands, and its 
energy eigenstates, magnetic Bloch states. A crucial difference between a Bloch state and 
a magnetic Bloch state lies in their translational properties. The Hamiltonian Hq is not 
invariant under lattice translation because Ao(r) cannot be a periodic function if the mean 
value of B is not zero. However, Hq can be made invariant under "magnetic" translation 
operators, which are the usual translation operators multiplied by a position dependent 
phase factor. @ 

We first give a brief review of the magnetic translation symmetry. In order to simplify 
the discussion, we assume the motion of electrons is confined in a plane (r = (x, y)) and the 
magnetic field is along the z direction. A magnetic Bloch state is the state that satisfies 

ffo^nk(r)=^„(k)^„k(r), (1.4) 

as well as 

Ti(i?i)vl/„k(r) = e''=^«^vl>„i,(r) 

T2(i?2)*nk(r) = e^'=^^^^„k(r), (1.5) 

where Ti and T2 are magnetic translation operators. Although Ti and T2 commute with the 
Hamiltonian by construction, they do not commute with each other unless there is an integer 
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number of flux quantum 0o enclosed by |Ri x R2I. Therefore, when the magnetic flux is a 
rational multiple p/q of the flux quantum 0o per unit cell of the lattice (plaquette), we must 
choose a "magnetic" unit cell containing q plaquettes in order that both ki and k2 be good 
quantum numbers, "^nk thus defined forms a complete set and satisfies the orthogonality 
condition 

{'^ n'\L'\^ nk) = Sn'nSk'k- (1-6) 

The domain of k is a magnetic Brillouin zone (MBZ), which is q times smaller than a usual 
Brillouin zone. Furthermore, because of the magnetic translation symmetry, the MBZ has 
exactly a g-fold degeneracy. We will call each repetition unit a "reduced" MBZ. 

One example of magnetic Bloch bands is the subbands split from Bloch bands due to a 
magnetic field. The number of subbands into which a band splits depends on the magnetic 
flux per plaquette in an intricate way.i If = p/g (in units of 0o), a Bloch band will split 
into q magnetic subbands. On the other hand, if the magnetic field is very strong, it is 
more appropriate to treat the lattice potential as a perturbation, then a Landau level will 
be broadened and split into p subbands. 

For a usual solid with a lattice constant a = 5A, the magnetic field has to be as large 
as 10^ Tesla in order for p/q to be of order unity. This is the reason the splitting was once 
considered impossible to observe. However, the field strength can be greatly reduced to 
a few Tesla if we use an artificial lattice with a much larger (say, 500A) lattice constant. 
Evidence for such splitting has appeared in recent transport measurements.! It is expected 
that more evidence will emerge in the future by using a very pure sample in a very low 
temperature environment. Under such circumstances, what is the dynamics for electrons in 
such magnetic Bloch bands? 

Using the magnetic Bloch states as an unperturbed basis, we found the following semi- 
classical dynamics in magnetic Bloch bands:0 

^=^-kxl7.(k), (1.7) 

and 

hk=-eE-erx6B, (1.8) 

where -E'n(k) consists of a band energy £^„(k) and a correction from the magnetic moment of 
the wave packet (this correction did not appear in Ref. 7). r2„(k) is the "Berry curvature", 
whose integral over an area bounded by a path C in k-space is the Berry phase r„(C)lE 
and are external fields added to the already present Bq field. These equations will be 
derived and explained in detail in Sec. II. 

Despite the similarities between Eqs. (1.1), (1.2) and Eqs. ( |1.7|) , (|1.8|) , there are several 
essential differences. See Table I for a comparison between this new semiclassical dynamics 
and the conventional one. The last item in Table I, about the quantization of orbits, will 
be explained in Sec. IV. We have to emphasize that the 6B in Eq. ( |1.8| ) is the field applied 
to the magnetic Bloch states; it is not the total magnetic field applied to the sample. This 
separation is particularly useful when B{x, y) is composed of a large constant part Bq and a 
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small non-uniform part SB{x,y). In this case, we can calculate the effect of Bq exactly and 
treat SB{x, y) classical perturbation. 

This paper is organized as follows: Sec. II is devoted to the derivation of Eqs. ( |1.7| ) 
and ( |1.8|) . Their use in calculating transport properties in a DC or AC electric field is 
demonstrated in Sec. III. The presence of 5B will lead to formation of cyclotron orbits in 
magnetic Bloch bands, similar to the formation of cyclotron orbits in usual Bloch bands. This 
is explained in Sec. IV. In Sec. V, we explore the connection between these cyclotron orbits 
and the Hofstadter spectrum. In Sec. VI, we estimate the energy levels in the Hofstadter 
spectrum by calculating the cyclotron energies in magnetic Bloch bands. Finally, this paper 
is summarized in Sec. VII. 



II. DERIVATION OF THE NEW SEMICLASSICAL DYNAMICS 



The method we use is to construct a wave packet out of \I'„k (hence it is already included 
the effect of i?o), and study its motion governed by the following Hamiltonian: 

H 2 



2m 



d 

-ih— + eAo(r) + e5A(r, t) 
or 



(2.1) 



where —d6A/dt = E, and V x 6 A = 5B. For simplicity we assume both E and are 
uniform; the derivation is still valid if they are slowly varying in space and/or time. 



A. Wave packet in a magnetic Bloch band 



Our derivation will be confined to one energy band by neglecting interband transitions; 
therefore, the band index n is henceforth dropped. Consider the following wave packet 
centered at which is formed from the superposition of magnetic Bloch states. 



\Wo) 



MBZ 



£kwik) |^(k)). 



(2.2) 



where w(k) is a function localized around kc (see Fig. 1 for an illustration). It has to be 
chosen such that 



and 



J d^k k |w(k)|2 = k. 



(Wo|r|W^o) = r^. 



(2.3) 
(2.4) 



By defining Mk(i") = e \Efk(i"), the mean position of Wo can be written as 



{Wo\r\Wo) = I d^k! j d^k w;*(k')w(k)(^(k') 



"9k^ 



ikr 



\u{k)) 



d^k' J d^k w*{k')w{k) 

d^k 



d 



d 



I— 5{k - k') + 5{k - k')(«(k)K— |«(k)) 



dk 



dk' 



cell 



w\k)^^w{k) + \w{k)\\u{k)\^^\u{k)U, 



(2. 
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FIGURES 




FIG. 1. Schematic plots for the motion of wave packets in r and k- space, 
where we have used the identitj;^ 

(«(kO|e^('^-'^')-M|^|«(k)) = 5(k-k')Kk)|z|^|«(k)),eii. (2.6) 

The subscript means that the spatial integration is restricted to a magnetic unit cell. By 
defining 

Ak)=2(M(k)||-|w(k))eell, (2.7) 

Eq. ( p.4|) can be written as 

d 



w*{k)t-^w{k) + \w{k)\^A{k) 



(2.8) 



B. Effective Lagrangian for a moving v^^ave packet 
The dynamics of a moving wave packet is governed by the following effective Lagrangian 

L(r„ k„ f-e, kc) = {W\th^^\W) - {W\H\W), (2.9) 

where IF is a wave packet centered at and k^ in the presence of external electromagnetic 
fields (kc is treated as a generalized coordinate here). We can always choose a gauge such 
that the vector potential SA is locally gauged away at a chosen point r = iv. At this 
particular point, the moving wave packet IF is the same as the Wo in Eq. (2.2). El The value 
of W near Tc can be approximated as 

IF(r) = e-''=/^^^(''-*)"VFo(r). (2.10) 
First, we evaluate the energy of this wave packet, which is (IF|if|IF) = (IFol-ff'lIFo), with 
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i^' = ^ I -ih-^ + e Ao(r) + e [5A(r, t) - 5A(r„ t)] | + V (r) 

^ + 7^ {[<5A(r, t) - 5A(r„ t)] ■ P + h.c.} . (2.11) 
2m 

P is the mechanical momentum operator corresponding to Hq. For simphcity, we choose the 
circular gauge for 6B, which gives 5A(r, t) = — Et + x r, and leads to 

H' Ho + —6B ■ L, (2.12) 
2m 

where L = (r — r^) x P is the mechanical angular momentum of the wave packet about its 
center of mass. Therefore, 

(W\H\W) ~ S(K) + —SB ■ (WolUWo). (2.13) 

2m 

The second term represents the energy correction due to magnetic moment of the wave 
packet. Notice that while a wave packet in an ordinary Bloch band does not rotate, a wave 
packet in a magnetic Bloch band usually does. 

For the first term on the right hand side of Eq. (2.9), we have 



{Wo\th^^\Wo) 



= {Wo\e6A{r„t) ■ r\Wo) + J (fk w*(k)i/i— w(k) 
= e(5A(re,t)-re + J d\ \w{k)\^h^p{k,t) 

~e(5A(r„t)-r, + ;i^7(k„t), (2.14) 

where w(k) = \w{k)\e~^"'^^'^\ Up to terms of total time derivative, which have no effect on 
the dynamics, the last line can be written as 

-eSA-h-nk,-^-f{K,t). (2.15) 



Using the condition in Eq. ( |2.8| ) and neglecting another term of total time derivative, we 
can write Eq. ( |2.15| ) as 

-e5A-r^ + hkc-rc + fi'kc-A{kc). (2.16) 

Combining Eqs. ( p.l3|) and ( p.l6| ), we have a final form for the effective Lagrangian 
(omitting subscript c) 

L(r, k, f, k) = -e5A(r, t)-r + nk-r + nA{k) ■ k - E{k), (2.17) 

where E{k) = S{k) + (e/2m)5B ■ L(k). Under a gauge transformation for Ao(r), the Berry 
potential -4,(k) will be changed by a term like (9koX(ko), and this will only change L by a 
total time derivative. This is also true if a different gauge is chosen for SB. Therefore, the 
dynamics is invariant under gauge transformation. 
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The dynamical equations in Eqs. (|1.7|) and (|1.8| ) can be obtained straightforwardly from 
this Lagrangian by using the Euler-Lagrange equation. The relation between r2(k) and 
A(k) is 

n{k) = V X A{k), (2.18) 
which can also be written as (z-component) 

This is the familiar Berry curvature in the study of the quantum Hall effect .0 

III. TRANSPORT IN MAGNETIC BLOCH BANDS 

A. Transport by an electric field 

The next step is to combine the semiclassical equations with Boltzmann equation to 
study the transport properties of magnetic Bloch electrons. The Boltzmann equation is 

r.^f^i^.'f^m , (3.1) 



dr (9k V dt 



coll 



where / = /(r,k) is a distribution function. The effect of impurities is included in the 
collision term {df/dt)^^^^. We use the relaxation time approximation to replace it by 
— (/ — /o)/T(k). For a random distribution of delta impurities voJ2i^{r — r^), the impurity 
scattering rate is 

riSp)-' = ^p{SF)n,vl (3.2) 

where rij is the area density of impurities. This rate is proportional to the density of states 
p{£f) at Fermi energy, which varies wildly with £p because the energy spectrum is discrete. 
However, since the following calculation is confined to only one band, r will be approximated 
by a constant. 

The equations of motion of an electron subject to a uniform electric field are 

where i^(k) is the reduced form of -E(k) in the absence of 5B. Substituting the expressions 
for r and k in Eq. (3.3) into Eq. (3.1), and setting / = /o on the left hand side of Eq. (3.1), 
we obtain 

where /o is the distribution function in equilibrium. Electric current is given by 
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It can be decomposed into three parts, + J"^ + J'", with the following definitions: 

J" = eV / — ^ I 1 Vb (vb ■ E) 







(27r)2 [ 
d^k / 




dfo] 


{2ny [ 


dS ) 



J'^ = er 77^ -1^ V v.^ . (3.6) 




In these expressions, Vb is the velocity dS/hdk due to energy dispersion, and v is the total 
velocity in Eq. (3.3). The meaning of these currents is explained below. 

First, is the Hall current. This is most evident considering a filled band with /q = 1. 
In this case, both and J'^ vanish, and only is nonzero. The integral of Q(k) over one 
magnetic Brillouin zone divided by 27r is always an integer, which is the topological Chern 
number discovered by Thouless et a/.0 Therefore we have 

= -C— E X z, (C G Z). (3.7) 

This formula represents the quantization of Hall current for a magnetic Bloch band. 

Second, 3^ is the diffusion current due to disorder scatterings. It can be put in the 
following form 

Vb(t)vb(0)-E, (3.8) 

where Vb(t) = e^*/'^Vb(0) is the current relaxed by scatterings after time t. At very low 
temperature, it can be simplified to 

3^ = e^g{SF)D-E, (3.9) 

where g^Sp) is the density of states at Fermi energy. D = c?^(vb(t)vb(0)) is the diffusion 
tensor, and the angular bracket ( ) means averaging over the Fermi surface.^ 

Third, J'^ is the current due to density gradient. It can be put in a form similar to 3^ in 
Eq. (3.9) at low temperature, but with two changes: (1) The velocity Vb in D is replaced by 
the total velocity v that includes the curvature term. (2) The driving force eE is replaced 
by dfi/dr. To evaluate this current, we need to know the explicit form of i7(k). This in 
general requires numerical calculation (see Sec. VI). 

We remark that, even though the derivation of Eq. (3.8) is based on a uniform electric 
field, its validity goes beyond that. It is actually a kinetic formulation, first proposed by 
Chambers, that is also valid in the presence of a magnetic field (for magnetic Bloch bands, 
the magnetic field is 55). Ill 
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B. Perturbation by an AC electric field 



The semiclassical method is much simpler to use than full quantum mechanical ap- 
proaches. This is most evident when the perturbation is changing in time. We illustrate this 
by considering a magnetic Bloch electron in an AC electric field. To simplify the discussion, 
we will neglect the effect of disorder and focus on the dynamics itself. 

Assuming that a uniform electric field along the x direction oscillates with a low frequency 
uj, we then have 

hk= -eEoe^^'x. (3.10) 
Considering a square lattice with the following energy spectrum 

S{k) = 2 [cos(A;ia) + cos(A;2a)] , (3.11) 
and substituting the solution k(i) into the r equation in (3.10), we have 

hx = 2a sin (- — — sm(ut) — kn^a 

cE 

hy = —2asm(kQya) —fl{kit)) cosi^ut) , (3-12) 

where (fcox, ^oy) is the initial value of k. It is not difficult to see that after many cycles of 
oscillation, there is a net drift along the x direction with average velocity 

{x) = — —sm{ko^a)Jo{z), (3.13) 

where Jp is the zeroth order Bessel function, and z = eEQa/hu is a ratio between two energy 
scales.cJ It can be seen that the original band transport velocity — 2a//isin(fcoxfl) is modified 
by Jo{z) because of the AC field. The electron is immobile along the x direction when 2 is a 
zero of the Bessel function. This resembles the collapse of the usual Bloch band in the AC 
Wannier-Stark ladder problem.ta 

Finally, we comment that the use of semiclassical equations is based on the assumption 
that impurities do not alter the band structure. Therefore, the electron dynamics between 
collisions can be nicely described by Eqs. (|1.?| ) and (|1.8|) . This assumption is no longer valid 
when q is large. In that case, disorder broadening tends to merge the subbands and wash 
out the fine structure. (This will be clearer after the discussion of hierarchical structure 
of the energy bands in Sec. V.) However, it was found that despite the energy spectrum 
has singular B dependence, the density of states appears as a continuous function of the 
magnetic field.0 Therefore, we divide total magnetic field B into Bq and SB, where Bq is 
related to the band structure undestroyed by disorder, and 6B is a small perturbation. In 
this case, the semiclassical dynamics in the magnetic Bloch band of Bq, driven by E and 
6B, will be employed in the Boltzmann equation. 
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IV. MAGNETIC PERTURBATION AND HYPERORBITS 



The usual Bloch electron will circulate around the Fermi surface along a constant en- 
ergy contour in the presence of a magnetic field. It is well-known that the quantization of 
cyclotron orbits leads to the famous de Haas- Van Alphen effect. In this section, we study a 
similar type of cyclotron motion in magnetic Bloch bands. It will be seen that this investi- 
gation yields very fruitful results. In particular, it offers a very simple explanation for the 
complex Hofstadter spectrum, which will be shown in Sec. V and Sec. VI. 



A. General properties of hyperorbits 

Combining the two equations in (3.3), we can eliminate r to obtain 

hk=-eZsB{k)^^xSB. (4.1) 
nok 

where ZsniM) = (1 + ^(k.)6Be/h)~^ is a curvature correction factor. This equation deter- 
mines the trajectories of magnetic Bloch electrons in k-space. It is not difficult to see that 
k moves along a constant energy contour of E(k) (which is slightly different from S{k)). In 
a classical picture, it is the drifting-center trajectory of the tighter cyclotron orbit formed 
from Bq. However, we have to emphasize that, the existence of hyperorbits is of quantum 
origin and cannot be explained classically. To differentiate them from the usual orbits of 
Bloch electrons, we will call them 'hyperorbits'. The hyperorbit in real space is derived 
from r = k X z{h/e6B), which is the k-orbit rotated by tt/2 and scaled by the factor h/e5B. 
It is also possible to define an effective cyclotron mass according to its frequency. However, 
this frequency will be very sensitive to its energy if the magnetic Bloch band is narrow, 
which is usually the case. 

There are several ways to verify the existence of hyperorbits. One way is through the 
measurement of magnetoresistance oscillation that originates from the quantization of hy- 
perorbits. This oscillation has a much shorter period than the usual de Haas-van Alphen 
oscillation because the effective magnetic field 6B is much smaller. The other way of veri- 
fying it is by using an electron focusing device to detect its real space orbit .0 This method 
has been used to map out the shape of a Fermi surface by measuring the shape of cyclotron 
orbits.Ei Another possible approach is to observe the ultrasonic absorption spectrum of the 
sample. The energy of an ultrasonic wave will be absorbed when it is in resonance with the 
hyperorbits. Similar method has been used to detect the existence of composite fermions in 
half-filled quantum Hall systems.0 



B. Quantization of hyperorbits 

In a previous paper, we have derived the quantization condition using Lagrangian for- 
mulation combined with path integral method.0 Here, it will be rederived using a slightly 
different approach. Substituting r = k x z{%/e5B) into the Lagrangian in Eq. (2.17), we 
will obtain an effective Lagrangian for the quasimomentum k. 
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fe2 

L(k, k) = ^(A;i^2 - hki) + hA-\i- E(k). (4.2) 
eoB 

It can be easily shown that Eq. (4.1) does follow from this L. The generalized momentum 
for coordinate k is equal to 

dL fi^ 

This leads to the following effective Hamiltonian 

i7(k,7r) = 7r-k-L(k,k) = ^(k). (4.4) 

Notice that because H does not depend on tt, the coordinate k will be a constant of motion 
and the dynamics is trivial. A Hamiltonian that gives correct dynamics will be given in 
Appendix A. Since it is not central to our derivation, we will not discuss it here. 

The quantization of hyperorbits is given by / tt -rfk = (m+7)/i, where m is a non-negative 
integer and 7 will be taken to be 1/2. It leads to area quantization in k-space 



- (f (k X dk) ■ z = 2tt (m + - 

2 J Cm \ 2 



where 

T{Cm) = <f A-dk (4.6) 

Cm 

is the Berry phase for orbit Cm- The orientation of Cm is chosen such that the sign of the 
area on the left hand side of Eq. (4.5) equals the sign of SB. 

The total number of hyperorbits in a MBZ is determined by requiring the area of the 
outer- most orbit be smaller than the area of a MBZ. Assume the flux before perturbation 
is BqO^ = p/q, then the number of hyperorbits in a MBZ is equal to |l/(g50) + where 
S(j) = 5Ba?e/h, and a is the Hall conductivity of the parent band. These hyperorbits are 
the lowest order approximation to the split energy subbands. They will be broadened by 
tunnelings between degenerate orbits. Since the MBZ is g-fold degenerate, the above number 
has to be divided by q to get the actual number of daughter bands, 

D = M?M±^. (4.7) 

Q 

This formula is essential in understanding the splitting pattern of the Hofstadter spectrum. 

The Hall conductivity for a subband can be calculated in the following way: In the 
presence of both E and 5B, the velocity of a magnetic Bloch electron consists of two parts,0 

• „ E X i , ^ ^ 

-kxz + ——. (4.8) 



e5B 5B 

The first term is the velocity of revolution, and the second term is the velocity of drifting 
along E X B direction. The current density for a filled subband is 
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J = - / d^kAk xz-e [ rf'k^-^. (4.9) 
J oB J on 

The first integral is zero for a closed orbit. Therefore, the Hall conductivity is obtained from 
the drifting term, which leads to (Jyx{oT simply a) = pe/6B, where p is the electron density 
per unit area. This will be used in the next section to determine the Hall conductivities for 
subbands in the Hofstadter spectrum. 



V. THE HOFSTADTER SPECTRUM 

A. Hierarchical structure of the spectrum 

The discussion in the preceding section presumes that we know 6B. But this is not 
apparent if the magnetic field B is homogeneous. In this case, we can still divide it into two 
parts, but where is the dividing point between Bo and 6B7 A natural way of dividing B (or 
(h) is to write it ClS db continuous fraction, 

0= , (5.1) 



fi + 



f2 + 



/3 + 



and truncate it according to the accuracy we need. The r-th order approximation of will 
be written as Pr /qr^ For example, if = 1/(2 + v^), we will have pi/qi = 1/3, P2/q2 = 
2/7, Ps/qs = 5/17- ■ ■ etc, which are truncations of 

(5.2) 



2 + V2 



2 + 



2 + --- 

The gr-'s satisfy the following recursion relation 

Qr+l = fr+lQr + Qr-l, (5-3) 

which relates the number of subbands at neighboring orders. There is also a relation between 
p's and g's, 

Pr+iq-r -Pr-Qr+l = ("I)''- (5-4) 

It follows that the extra magnetic flux between the r-th order and the (r + l)-th order 
truncation is 

50. = ^ - ^ = tlL, (5.5) 

Qr+l Qr Qr+lQr 

Since the sign of 50. alternates from one order to the next, the direction of 6B also 
alternates.ii 
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Notice that for a chosen fraction Pr/<iri the size of a MBZ is fixed. Without perturba- 
tion from the part that is truncated away, a wave packet will move on a straight line. The 
trajectory is curved because of 50^- Higher the order approximation we use, smaller the 50 
gets and larger the radius of the hyperorbit becomes. In the ideal case, without any compli- 
cations due to disorder, thermal broadening ■ ■ ■ etc, we expect there will be a hierarchical 
structure of hyperorbits due to different orders of approximation. This structure finds its 
correspondence in the hierarchical structure of the Hofstadter spectrum.ii (Fig. 2) 
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FIG. 2. Hofstadter spectrum for a square lattice generated from a tight-binding model (g < 60). 
Only the < 1/2 part is shown since the spectrum is symmetric with respect to = 1/2. 



B. Distribution of Hall conductivities and splitting of energy bands 

A magnetic Bloch band carries quantized Hall current, and this current will redistribute 
among daughter bands in such a way that the total Hall current for subbands equals the orig- 
inal current.^ We will call this the "sum rule". The current distribution among subbands, 
which is also quantized in each subband, was obtained by Thouless et a/.0 In their famous 
paper, they found that the subband Hah conductivities are the integer-valued solutions of 
the Diophantine equation. Here we show that semiclassical dynamics offers an alternative 
and very heuristic solution to this problem. The Hall conductivities calculated will be used 
in Eq. (4.7) to calculate the number of magnetic subbands after splitting. 
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The general expression for the Hall conductivity of a "closed" subband is o" = pe/5B (see 
Eq. (4.9) and below). Therefore, a can be determined if p and 6B are known. Consider a 
subband at the r-th order of splitting. Since all subbands at this level share the same number 
of states, each subband will have Pr = Po/Qr, where po is the density of states for the original 
Bloch band. The perturbation field for a subband at this level is 5Br — hS(f)r-i/{ea^). 
Therefore, we have 

af°'^ = |gl = (-i)-ig^_,, (in units of e^h). (5.6) 

(Since oi — 1 for a closed subband at the first level, qq will be set to one .) We have to 
emphasize that, Eq. (5.6) is valid for every closed subbands at the r-th order. Combining 
Eq. (4.7) with (5.6), we can determine the number of daughter bands being split from a r-th 
order parent band, which is 



close 



lose I 



Ui. (5.7) 



FIG. 3. Schematic plots of the constant energy contours in the reduced MBZ of a square lattice 
and a triangular lattice. Dashed lines are the open orbits. 

The Hall conductivity for an open subband is more difficult to obtain. It requires the 
knowledge of the exact k-trajectory to figure out the first integral in Eq. (4.9). However, for 
a square or a triangular lattice, there is an easier way of calculating it. This is so because, 
there is only one open orbit in every MBZ for either lattice (see Fig. 3). Therefore, there 
is only one open daughter band for every parent band. Its Hall conductivity can be figured 
out by using the "sum rule" : 

(^parent = '^daughter ■ (5-8) 



For example, at the first order, we have a°^^^ — 
since it is a Bloch band at the very beginning, 
daughter band at the r-th order is 



<Jo - ifi - l)<7f°^^ = -(/i - 1), where (7o = 
For r > 2, the Hall conductivity of an open 



—Open ^close / / -\\ ^c 
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= (-irV2 - (/. - i)(-ir Vi 



(5.9) 



where we have assumed that its parent is a 'closed' band. Using Eq. (4.7), we see that this 
open band (now as a parent) will split into 

Pr" = /r+l + l (5.10) 

subbands under perturbation. It can be shown that the same result as Eq. (5.9) is obtained 
if its parent is an 'open' band (with /r + 1 daughters). This checks the consistency of this 
calculation. It is clear that the extra splitting of one subband from each open parent band 
(there are qr-i of them) accounts for the extra Qr-i in the recursion relation Eq. ( ^.3|) 

We give one example to demonstrate the use of these rules. Consider a square lattice 
with = 1/(2 + v^)- Because af°^^ = 1, and 0"°'''''^ = —2, the distribution of cr's for the 
three subbands at the first order is 



0-1 = (1,-2,1). (5.11) 

where we have put 0"°^'^" in the middle since for a square lattice the open subband is located 
at the center of a parent band (Fig. 3). 

These three bands will be spht into seven subbands due to the extra flux S(f>i =^2/^2 — 
Pi/Qi = 2/7 — 1/3. Since /2 = 2, the pattern of splitting will be, according to Eqs. (5.7) and 
(5.10), 

X>i = (2,3,2). (5.12) 

Furthermore, since af°^'^ = —3 and a2^'^^ = 4 according to Eqs. (5.6) and (5.9), the Hall 
conductivity distribution is0 

0-2 = (-3,4,-3,4,-3,4,-3). (5.13) 

Consequently, we have 

X>2 = (2,3,2,3,2,3,2). (5.14) 

The Hall conductivities we just obtained are the same as those derived from the Diophantine 
equation. Actual pattern of splitting is shown in Fig. 4 for comparison. The distribution in 
Fig. 4 for the left (or right) five subbands in Eq. (5.14) appears to be (2, 1, 2), instead of (2, 3) 
(or (3, 2)). However, closer examination reveals that the left (right) three subbands actually 
come from the same parent. In fact, slight asymmetry in the distribution is inevitable 
because when is changed by a small amount, an electron state cannot suddenly jump out 
of the band edge to the middle of a gap. 

Eqs. (5.6)-(5.10) also apply to a triangular lattice. The only difference is that (Xopen no 
longer locates at the center of a parent band. Given the same = 1/(2 + -\/2), we now have 

Ti = (l,l,-2) 
X>i = (2,2,3) 

0-2 = (4,-3,-3,4,4,-3,-3) 
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X>2 = (3,2,2,3,3,2,2). (5.15) 

This again conforms with the actual spectrum and the solutions of the Diophantine equation 
with subsidiary constraints suitable for a triangular lattice.E^ 



2/7 



FIG. 4. Pattern of energy splitting for (j) in Eq. (5.2). Different orders of approximation give 3, 
7, and 17 subbands respectively. 

VI. CALCULATIONS OF ENERGY SPECTRUM, CURVATURE, MAGNETIC 
MOMENT, AND CYCLOTRON ENERGY 

In this section, we give detailed calculations of Sn(X), ^n(X), and Ln(k). They are used 
in calculating the cyclotron energies using the quantization formula Eq. ([4.5|) . The cyclotron 
energies will be used to estimate the subband energies in the Hofstadter spectrum. In doing 
so, we not only presume the one-band approximation (on which Fig. 2 is based), but also 
neglect the inter-orbit transitions that broaden (and may slightly shift) the energy levels. 
The latter approximation leads to negligible error if the bandwidths under consideration are 
very small. Similar calculations have been done by Wilkinson, and his results have been very 
successful.0 However, we believe that the approach proposed here is conceptually simpler, 
and is easier to generalize to other types of lattices. 



A. Calculation of energy spectrum for parent bands 

The following calculation is based on the tight-binding model. We will only give a very 
brief explanation of this approach. For more details, we request readers to refer to Ref. (29). 
In the tight-binding approximation, a Bloch state is expanded as (for (f) = p/q) 

M/(k)=^aKk)^Kk), (6.1) 

1=1 

where k is restricted to a reduced MBZ. The basis ipiik) is defined to be ipi^i + 27r0/, ^2), 
where V'(k) is a Bloch state before the magnetic perturbation. A tight-binding Hamiltonian 
in the absence of 6B, when being expressed on the basis of is a g x g matrix (the lattice 
constant a is set to one.) 
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/ cos(fci + 27r/g) 



I' I 



cos(A;i + 47r/g) e 





ik2 



\ 



-ik2 











□ «fc2 



cos(A;i + (g-l)^^ 

g«fc2 






giA;2 

cos(A;i) y 



(6.2) 



£^„(k) and a['(k) are nothing but the eigenvalues and eigenvectors of this matrix. For exam- 
ple, if p/q = 1/3, then a straightforward calculation shows that the £„(k)'s are solutions of 
the following characteristic equation. 



£^ + 6£ = 2 [cos(3A;i) + cos(3A;2)] . 



(6.3) 



There are three roots for each k, and variation of k over the MBZ leads to the three energy 
bands for = 1/3 (see Fig. 2). It is not difficult to see that the band edges are located at 
(from high to low) Si{g) , Si{0) , 82(0) , S2{g) , S^ig) , and £:3(0), where g = (7r/3,7r/3). 



B. Calculation of Berry curvature 



To calculate the Berry curvature in Eq. ( |2.19| ), we need to know the eigenvectors of Hin. 
Before doing that, we will try to rewrite Eq. ( p.l9 ) in a form suitable for the tight-binding 
calculation. First, we insert a complete state J2n' \un'){un'\ inside the dot products that 
appear in Eq. ( |2.19|) . Since we are using the one-band approximation, n' only runs through 
subbands in the same parent band, and 



fi„(k) = z£ ' 



n'=l 



{^^\Un'){Un' 
OKi 



dUn 
^dko 



) — c.c. 



(6.4) 



where we have dropped a term with n' = n since {un\d/dk\un) is purely imaginary and does 
not contribute to the curvature. With the help of the identity 



/ I ^ I \ (M„'||k|Mn) 



'(9k' 



(6.5) 



where H = e ^'^ ''He (this and the following iJ's are the unperturbed Hamiltonian in 
Eq. (1.3); the subscript '0' is dropped for brevity), we can rewrite Eq. (6.4) in the form 



fi„(k)=z^ ' 



dH I 



— C.C. 



Expanding m„ by ipiklr), which is defined to be e~'^^'^ipik{r), we have 



(6.6) 



{Un' 



OH. 
'dk 



\u 



~ .dH ~ 



dk 
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+ {£n-£n')Y:atanM^)- (6.7) 

The inner product {ipi'\dipi/dk.) is zero for a Bloch state in the absence of a magnetic field. 
Therefore, 

{un'\^\un) = l^ai O/^j^- (6.8) 

Beyond this stage, the calculation is straightforward since we only need to calculate af from 
Eq. ( |6.2| ) and combine Eqs. ( |6.6| ) and (6.8) to obtain f2„(k). 

Again we choose a simple fraction p/q = 1/3 and calculate the Berry curvature distri- 
butions for the three magnetic subbands. The result is shown in Fig. 5, in which the range 
of k vector is one reduced MBZ (the basic unit of repetition). Note that the curvature 
tends to concentrate on four inner band edges because the electron states near inner gaps 
are changed the most from the original Bloch states that have zero Berry curvature. The 
curvatures from the three bands cancel locally, that is 

3 

^a„(k) = 0, Vk. (6.9) 

n=l 

This is in general true for any q and can be easily proved from Eq. ( |6.6| ). Finally, integration 
of flniX) over a MBZ divided by 27r gives us integers (1, —2, 1). These are indeed the Chern 
numbers we expected (see Eq. (5.11)). 
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FIG. 5. Distributions of Berry curvature ^in(k) (in units of /h). r2i(k) is equal to J^3(k) 
shifted by (7r/3,7r/3) 



C. Calculation of magnetic moment 

The self-rotating angular momentum of the wave packet (L)„k in Eq. (2.13) can be 
written in a from that is more tractable for calculation: 



i:„(k) 



(6.10) 



Derivation of this formula is given in Appendix B. Eq. (|6lO| ) can also be rewritten as 



.m 



C.C. 



(6.11) 



Derivation of Eq. ( p.ll| ) is very similar to the derivation of Eq. (|6.6|) . The only change is 
that the extra factor of H — Sn in the numerator cancels a Sn' — £n in the denominator. 

I/„(k) can be readily calculated by combining Eq. ( |6.8| ) with Eq. ( |6.11| ). The result is 
shown in Fig. 6 (again for = 1/3). Similar to the distribution of f2„(k), L„(k) also has 
peaks near inner band edges. The total magnetizations from all three bands cancel each 
other. However, unlike the curvature, they do not cancel locally at each k point. 
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We remark that this magnetization energy first appeared in a paper by Kohni. Their ob- 
jective was to study the effective Hamiltonian for Bloch electrons in a weak electromagnetic 
field. This term is in general zero for Bloch bands, but can be nonzero in the presence of 
spin-orbit interaction. In the latter case, it contributes an extra g-factor to Bloch electrons.ii 
An expression that is the same as the right hand side of Eq. ( |6.10| ) has also been obtained 
by Rammal and Bellissard.0 Without the m/h factor (and apart from a factor of two), it is 
called the Rammal- Wilkinson form. 






FIG. 6. Distributions of angular momentum L„(k) (in units of 2m /fi). 
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D. Calculation of cyclotron energy 



After obtaining £^„(k), r2„(k), and L„(k), we can determine the cyclotron energies accord- 
ing to the quantization formula Eq. ( [4 .51) . In the following example, we add 5(f) = —1/201 
to (p = 1/3- This gives (p' = 22/67 = 1/(3 + 1/22). According to the simple rules derived 
in previous section, the original magnetic bands are expected to split to 22, 23, and 22 
subbands respectively. In Table II, we compare the exact spectrum with the quantized cy- 
clotron energies. Only part of the subbands from the parent band in the middle are shown. 
We show two numbers (for two band edges) in the first row when the bandwidth for -EHofst 
is larger than 10~^. The i^cycio's in the second row are obtained by fine-tuning the path 
in Eq. ( [4.5|) , with uncertainty on the order of 10~^. It can be seen that the match between 
ii^cycio and -Enofst is quite satisfying. We have done calculations for subbands from other 
parent bands, and they also show similar accuracy. 

Notice that the energy levels from fractions like = 1// are broadened cyclotron levels 
in an ordinary Bloch band. They do not split from a magnetic parent band. In this case 
f2„(k) and L„(k) are zero, and Eq. ( [4.5|) reduces to the usual Onsager quantization formula. 
Numerical result based on this simplified formula for the cyclotron energies also agrees very 
well with the positions of subbands in Fig. 2. 

VII. SUMMARY 

Electron states in a lattice subject to a homogeneous magnetic field satisfy magnetic 
translation symmetry and have band-like energy spectrum similar to the usual Bloch band. 
However, to our knowledge, the semiclassical dynamics of magnetic Bloch electrons has never 
been an explicit subject. One reason is that the observation of band-splitting remains an 
experimental challenge to date; the other reason might be due to the fact that magnetic 
Bloch bands, unlike Bloch bands, can be changed easily by varying an external magnetic 
field. However, our study has shown that the inquiry of magnetic Bloch bands can be very 
rewarding in itself. Major findings in this paper are summarized below: 

The Berry curvature of magnetic bands plays a crucial role in the dynamics. It gives 
electrons an extra velocity in the direction of E x B, and this term directly relates to the 
quantization of Hall conductivity. This semiclassical dynamics, combined with the Boltz- 
mann equation, is used to study electron transport in a DC or AC electric field. 

In the presence of 6B, the energy dispersion -E(k) is shifted from the usual band energy 
because of the non-zero magnetic moment. Similar to usual Bloch electrons, magnetic Bloch 
electrons execute cyclotron motion on the constant energy surface of -E(k). However, the 
quantization condition for cyclotron orbits has to be modified from the usual Onsager con- 
dition because of the Berry phase. Based on this modified formula, we obtain a simple rule 
that calculates the number of daughter bands for every parent band in the Hofstadter spec- 
trum. Furthermore, a fairly heuristic explanation for the distribution of Hall conductivities 
is given using the picture of cyclotron orbit drifting. 

These quantized orbits are closely related to the energy bands in the Hofstadter spectrum. 
We give detailed numerical calculations for f2n(k) and Ln(k) based on the tight-binding 
model for the case of = 1/3. They are used in the calculation of cyclotron energies 
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using the quantization condition, and the result is in very good agreement with the actual 
spectrum. This shows that the complex pattern of the Hofstadter spectrum is nothing more 
than the broadened cyclotron energy spectrum in magnetic Bloch bands. 
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APPENDIX A: 

Note that there is no k dependence in tt (because there is no kinetic energy term in 
L(k, k)). In this case, it is more pertinent to treat Eq. (4.3) as a constraint on the variables 
k and tt, 

6>(k,7r) = 7r + k X z - hA(k) 0. (Al) 

2eoB 

Strictly speaking, this constraint cannot be used before obtaining the equations of motion, 
therefore we use ~ to distinguish it from a real identity.S A general Hamiltonian for a 
system with constraints is given by, 

H*{k, tt) = H{k, 7z) + X-0 = E{k) + A • 6>, (A2) 

where A = (Ai, A2) are arbitrary functions of k and tt. The dynamical equations for this 
new Hamiltonian are 

dH* dE ^ ^ ^ ^ dAi 

where we have discarded a term dX/dk -6 = 0. According to Eq. (4.3), we should also have 
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-kxz + hY^k,—. (A4) 



2e5B Y dki 

Equating (A3) to (A4), and replacing A by k, we will get Eq. (4.1). 



APPENDIX B: 



We will rewrite the angular momentum of a wave packet in Eq. (2.13) in terms of magnetic 
Bloch functions. By defining w{k) = e^^'^^wlk), we have 
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L„(ke) = J d'k'J d^k ^*(k')w;(k)(^„(k')|(r-r,) x P|^„(k)) 

d^k' [ d^kw*{k')w{k){un{k')\e'^''-'''>^''''^\r~r,) x P(k)|M„(k)), (Bl) 



(B2) 



where P is the momentum operator on the |-u„) basis. Since 

p(k)|M„(k)) = |M„,(k))(M„,(k)|p(k)|M„(k)), 



and 



{un{k')\e 



'^U«(k-k')-{r-r,)/ 



rc)|M„/(k)) 



d 



i6n 



•"9k 



-6{k-k')-z6{k-k'){^M 



we have 



Ln(kc) 



z d'k 



*(k) X (P)„ 



,du. 



d'k\w{k)\'{^\xP{k)\Ur.). 

ok 



Because P(k) = {'m/h)dH /dk, the integrand of the second term can be written as 
m.dun:dH 

m d .dun.fj. > m dun.fj.dun. , ^ 



h dki dk2 h dki dk2 

^(^|x(^„-if)|^) + (^K)x(P) 



(B3) 



(B4) 



(B5) 



Therefore we have 



-i jd-'k |^;(k)p(^K) X (P), 



i / d'^k 



d. 



?2;(k) X (P)„. 



(B6) 



The last two terms cancel because of Eq. ( p.8|) , and this leads to Eq. ( |6.10| ) Notice that 
this result is independent of the way a wave packet is constructed since there is no w{k)- 
dependence in the new expression. 
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TABLES 



TABLE L 



Bloch band 



Magnetic Bloch band {{p/q)^o per plaquette) 



Unperturbed Hamiltonian = V ^(r) Ho = ^ [-in-§^ + eAoir)y + V (r) 



Translation operators 

Number of plaquettes 
per unit cell 

Range of k vector 



Perturbing fields 
Velocity of electron 

Dynamics for k 

Quantization condition 
for cyclotron orbits 



r(R) = 



1 plaquette 



One Brillouin zone 



E,B 



r = d6nOsi)/hdlc 



-eE — er X B 



q plaquettes 



One magnetic Brillouin zone 
(One Brillouin zone divided by q) 

E,(5B 

f = dEn{k)/hdk - k X n„(k), 
En{k) = f-^g(k) + {e/2m)6B ■ L„(k) 



hk 



-eE — er X 5B 



Area(C^) = 27r (m + i) eB/h Area{Cm) = 27r (m + | - ^^^) e6B/h 



TABLE IL Hofstadter spectrum (for (p = 22/67) and the cyclotron energies calculated from 
Eq. (4.5). Only the top ten subbands for the middle parent band are shown. The last column is 



the suV)} 


)aud ck)S(^sl lo Wic 




■)aiid edge 


^2(0) = 


0.7321. 










-E'Hofst 


0.0618 


0.1067 


0.1566 


0.2124 


0.2747 


0.3443 


0.4221 


0.5098 


0.6098 


0.7266 




0.0678 


0.1085 


0.1570 


0.2125 














-^cyclo 


0.0632 


0.1063 


0.1558 


0.2115 


0.2738 


0.3435 


0.4212 


0.5086 


0.6082 


0.7240 
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